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ABSTRACT 

Context. Convection in stars excites resonant acoustic waves which depend on the sound speed inside the star, which in turn depends on 

properties of the stellar interior. Therefore, asteroseismology is an unrivaled method to probe the internal structure of a star. 

Aims. We made a seismic study of the metal-poor subgiant star v Indi with the goal of constraining its interior structure. 

Methods. Our study is based on a time series of 1201 radial velocity measurements spread over 14 nights obtained from two sites, Siding 

Spring Observatory in Australia and ESO La Silla Observatory in Chile. 

Results. The power spectrum of the high precision velocity time series clearly presents several identifiable peaks between 200 and 500 /iHz 
showing regularity with a large and small spacing of Av = 25. 14 ± 0.09 /jHz and 6v 02 = 2.96 ± 0.22 yuHz at 330 /iHz. Thirteen individual modes 
have been identified with amplitudes in the range 53 to 173cms~'. The mode damping time is estimated to be about 16 days (1-cr range 
between 9 and 50 days), substantially longer than in other stars like the Sun, the a Cen system or the giant f Hya. 
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1. Introduction 

The analysis of the oscillation spectrum provides an unrivaled 
method to probe the stellar internal structure. The frequencies 
of these oscillations depend on the sound speed inside the star, 
which in turn depends on density, temperature, gas motion and 
other properties of the stellar interior. High-precision spectro- 
graphs have led in recent years to a rapidly growing list of solar- 
like oscillation detections. 
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The recent achievements are reviewed by Bedding & 
Kjeldsen (2006, 2007 ). All of the previous detections of oscil- 
lations were made on stars with metallicities close to solar. We 
recently reported observations of the metal-poor subgiant star 
v Indi made in both Chile with the spectrograph CORALIE and 
in Australia with UCLES (Bedding et al. l2006l hereafter Paper 
I). These observations covered a 14 nights time span and had 
the benefit of two-site coverage. 

In Paper I, we described the radial velocity time series 
and calculated the corresponding power spectrum weighted by 
the measurement uncertainties. We reached a noise floor of 
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Fig. 1. Power spectrum of the radial velocity measurements of v Indi. The inset shows the spectral window. 
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14.9 cm s 1 in the combined amplitude spectrum, which is quite 
high compared to other studies with the same instruments (see 
e.g. Bouchy & Carrier [20021 and Bedding et al. 120041 . This 
difference can be explained by the extremely low metallicity 
(about 3% of solar) and the faint magnitude (my = 5.28) of 
v Indi. Oscillation modes are present around 0.32 mHz with a 
maximal amplitude per mode of 95 cms -1 . A large frequency 
separation was determined by autocorrelation and Bayesian 
methods to be 24.25 + 0.25 //Hz. 

In this paper, we will extend the analysis by extracting indi- 
vidual oscillation frequencies, by measuring the large and small 
separations, and by estimating the mode lifetimes. 

2. Frequency analysis 

The power spectrum of the time series, shown in Fig. [1] ex- 
hibits a series of peaks between 200 and 550//Hz. In solar-like 
stars, p-mode oscillations are expected to produce a character- 
istic comb-like structure in the power spectrum with mode fre- 
quencies v n ,( reasonably well approximated by the asymptotic 
relation (Tassoul [l980l l: 



Vnjt * Av(n + -+€)-£(£ +l)D . 



(1) 



Here Dq (which equals ^6vq2 if the asymptotic relation holds 
exactly) and e are sensitive to the sound speed near the core and 
in the surface layers, respectively. The quantum numbers n and 
I correspond to the radial order and the angular degree of the 
modes, and Av and Svo2 are the large and small spacings. To 
fit to this relation, the large separation was measured in Paper I 
by autocorrelation and Bayesian methods and was found to be 
24.25 +0.25 //Hz. 

Note that a subgiant such as v Indi is expected to show 
substantial deviations from the regular comb-like structure de- 
scribed above. This is because some mode frequencies, except 
for I — 0, could be shifted from their usual regular spacing by 
avoided crossings with gravity modes in the stellar core (also 
called 'mode bumping') (see e.g. Christensen-Dalsgaard et al. 



|1995l and Fernandes & Monteiro 2003). We must keep in mind 
the possibility of these mixed modes when attempting to iden- 
tify oscillation modes in the power spectrum. 



Frequency modulo 24-25 /iHz Frequency modulo 25-4 /-iHz 

Fig. 2. Echelle diagrams for two different values of the large 
separation. The modes with a SNR greater than 4 are indicated 
with a square, and the ones with a SNR greater than 5 are in- 
dicated in boldface. Left: Echelle diagram where 24.25 //Hz is 
used as a value for the large separation (see Paper I). Right: 
Echelle diagram with a large separation of 25.4//Hz. The lines 
are just meant as a "guide". 



The next step is to measure the frequencies of the strongest 
peaks in the power spectrum. The frequencies were extracted 
using an iterative algorithm which identifies the highest peak 
between 150 and 600 //Hz and subtracts it from the time se- 
ries (see e.g. Carrier et al. [2003 and Kjeldsen et al. 2005]). In 
the present analysis we are only marginally affected by the 
window sidebands compared to mono-site observations anal- 
ysis. The highest amplitude modes are therefore expected to 
be located at their correct position and they do not need to be 
shifted by the daily aliases (Kjeldsen et al. 2003, Carrier et al. 
2005a 2005b). In total, 27 frequencies were extracted having 
an amplitude 3 times greater than the noise. These amplitudes 
are measured from the best-fitting sinusoid, whereas the noise 
is computed as the mean in the amplitude spectrum, which is 
the square root of the power spectrum, at high frequencies. 
Note that the noise level at high frequencies (600- 1000 //Hz) 
has a value of 14.9 cm s _1 (Paper I) and is assumed constant on 
the whole range 150-600 //Hz. 
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All frequencies with a signal-to-noise ratio (SNR) above 3 
were then placed in an echelle diagram, without shifting any of 
them, constructed with the value for the large separation found 
in Paper I (see Fig.[2]left). This diagram does not seem to show 
a classical asymptotic structure, since only one vertical ridge 
can be seen (dotted line at 22 yuHz). A possible way to explain 
such a graphic is to associate the ridge at 22yuHz to t — 2 
modes and the quasi-line near 27 //Hz to I = modes. This 
identification is not very convincing and would imply a small 
separation of 5 //Hz, which is too large for this kind of star. 

Another way to understand this diagram is to increase 
slightly the large spacing and allow curvature along the fre- 
quency axis (see Fig. [2] right). Indeed, large separation deter- 
mination methods based on the asymptotic relation (as the au- 
tocorrelation presented in Paper I), do not take into account any 
curvature. With curvature, there is no correct value for the large 
spacing, it is a function of frequency. With a large separation 
of 25.4 //Hz, as in the right panel, we can identify ridges corre- 
sponding to t — and 1-2. These are marked with dashed 
lines. 

With the ridges for I — and 2 marked as shown, we then 
used the asymptotic relation to calculate the expected positions 
for I = 1 and 3. These are shown by dotted lines in the right 
panel of Fig. [2] We can see that two strong peaks lie exactly 
on the line for £ = 3, giving us some confidence that we have 
found the correct solution. We also see an absence of peaks 
on the £ = 1 line and a few strong peaks that do not lie on 
any of the lines. This is consistent with the presence of mixed 
t = 1 modes that are shifted by avoided crossings, as described 
above. 

Based on Fig. [2j several modes were identified in v Indi. 
The frequency, amplitude, SNR of these modes are given in 
TableQ]with also the identified degree. Since the amplitude of 
the modes in an iterative fitting algorithm depends slightly on 
the selection order, a least squares fit was applied at the end to 
determine all amplitudes simultaneously. We retained all peaks 
with a SNR greater than 4. We also added two weaker peaks 
that were in good agreement with the t = 2 ridge (see TableHJ. 
It is expected in evolved stars that I = 1 modes are the first 
to be mixed modes and do not follow the asymptotic relation. 
Thus, all large peaks which can not be determined as t — 0, 2 or 
3 are listed as I = 1 modes. I — 3 modes are also less reliable. 
Finally, all the selected modes are shown in the echelle diagram 
on Fig. [3] and listed in Table [2] 

Based on Fig. [3] and on the frequencies given in Table [2] 
the large and small separation are deduced taking into account 
only the t — and 2 modes, as { = 1 modes do not follow the 
asymptotic relation and I = 3 ones do not have a secure identifi- 
cation. The large separation has a value of 25. 14 + 0. 14 //Hz at 
0.33 mHz which corresponds to the weighted average of indi- 
vidual values given for each degree (see Table|2]i. The variation 
of the large separation versus frequency is presented in Fig. |4] 

For the small separation, there are no pairs of adjacent 
I — 0, 2 modes, and so missing I — modes were linearly inter- 
polated. We find a value of 6v QZ = 3.0 + 0.3 //Hz at 0.33 mHz. It 
is not possible to calculate a reliable value for 5vo since C=\ 
modes do not at all follow the asymptotic relation, and since 
the identification of I — 3 modes could be in error. In fact, one 



Table 1. Amplitude and frequency for individual oscillation 
modes and their identification. The mean noise in the ampli- 
tude spectrum is 11. 9 cms -1 . 
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Fig. 3. Echelle diagram for the frequencies given in Table [T] 
and|2] The lines are just meant as a "guide". Squares, crosses, 
rhombuses and triangles represent the £-Q, 1, 2 and 3 modes 
respectively. 

could also identify the £-3 as £=l modes. This fact should 
be remembered when comparing the observed and theoretical 
frequencies. 

3. Mode lifetimes 

We have used two methods to estimate the mode lifetimes from 
our data. The first uses the distribution of mode amplitudes and 
the second uses the scatter of the frequencies about the ridges 
in the echelle diagram. In both cases, we calibrated the results 
using simulations. 
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Table 2. Frequency (in //Hz) and mode identification for in- 
dividual oscillation modes. The errors on the frequencies are 
given in brackets (10 _1 //Hz) according simulations described 
in Section [3731 The identification of frequencies in brackets is 
less secure. 
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Av 


25.14 
± 0.22 




25.14 
±0.18 





25. 14 ± 0. 14 yuHz at 0.33 mHz 



3. 



Frequency (juHz) 

Fig. 4. Large spacing Av versus frequency for p-modes of de- 
gree I = (□) and I — 2 (o). Open symbols correspond to large 
spacing averages taken between non-successive modes. 



3.1. Mode lifetimes by amplitude comparison 

The idea is to compare the amplitudes of peaks in the power 
spectrum with simulations that have the same amplitude enve- 
lope as the observations and a range of mode lifetimes. The 
simulator for generating solar-like oscillation time series is 
similar to that described by De Ridder et al. (2006). For the 
input amplitudes, we used the smoothed amplitude spectrum 
show in Fig. 7 of Paper I. We ran ten simulations for each of 
six mode lifetimes (1, 2, 5, 10, 20 and 50 days). Two examples, 
with mode lifetimes of 2d and lOd, are shown in Fig. [5] We 
see a significant difference between them, which can be used 
to estimate the mode lifetime. 

In order to do this, we extracted the 7 highest peaks from 
each simulated power spectrum (using iterative fitting) and 
compared their amplitudes with those extracted from the ob- 
served series. As a simple measure of the degree of agreement, 
we adopted the following method. For a given simulation, we 
compared the highest peak with the highest peak in the obser- 
vation. We then compared the second highest simulated peak 




400 600 
Frequency [/^Hz] 

Fig. 5. Examples of two simulations of the v Indi time series us- 
ing mode lifetimes of 2 days (top) and 10 day (bottom). Higher 
peaks are in general found for the longer lifetime. 

with the second highest observed peak, and so on. The fraction 
of peaks in the simulation that were above the observed ones, 
in this pairwise comparison, gave the measure of agreement, 
with 50% indicating the best agreement. 

The results are shown in Fig. [6] which shows the fraction of 
simulated peaks that were above the observations. We can read 
off the most likely value of the mode lifetime, corresponding to 
a fraction of 0.5, to be 17 days. However, this measurement is 
not very precise. The dashed curve is the cumulative distribu- 
tion of a normal distribution, and we see that the one-cr range 
for the mode lifetime is 2 to 141 days, or in the logarithmic 
units: 

log 10 (T/ d) = 1.23 + 0.92. (2) 



3.2. Mode lifetime by frequency scatter 

The other way to extract information on mode lifetimes is to 
examine the scatter of measured frequencies about the ridges 
in the echelle diagram. For these observations, this is only pos- 
sible for the five £ — modes. If the modes were coherent 
(i.e. pure sinusoids with infinite lifetimes), we would be able to 
measure them with a typical accuracy of 0. 1 /iHz. The observed 
scatter is much greater than this, which we interpret, at least for 
i — modes, as reflecting the finite lifetimes of the oscillation 
modes. 

As is well established for the Sun, the power spectrum of 
a stochastically excited oscillation that is observed for long 
enough will display a series of closely spaced peaks under a 
Lorentzian profile, the width of which indicates the mode life- 
time (e.g. Toutain & Fr6hlich [1992| ). If the observations are not 
long enough to resolve the Lorentzian profile, then the effect 
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Fig. 6. Calibration of the simulations to determine mode life- 
time from the peak amplitudes. It shows the fraction of modes 
in the simulations that have amplitudes above the observed 
mode amplitudes. The dashed curve is a normally distributed 
function that presents the best fit. Vertical dashed lines show 
the resulting mode lifetime and the values at ±l<x. 
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Fig. 7. Relation between scatter per mode and mode lifetime in 
days. The curve is a fit to points that came out of 700 individual 
simulated time series data. The mode lifetime corresponding to 
the best fitted frequency scatter is shown together with values 
for ±\cr. 



mode could be calculated, which have a mean of 0.36yuHz + 
0. 15 //Hz, for an average observed amplitude of 95 cm s . 

In order to calibrate the relation between frequency scat- 
ter and mode lifetime, we ran 700 simulations using the same 
sampling and statistical weighting as the real time series. The 
results are shown in Fig. [7] and the best fit corresponds to a 
mode lifetime of 16.2 days (between 8.7 and 60.8 days at 1 cr). 
Expressed in logarithmic units, we get: 

log 10 CrA0 = 1.21^f 7 (3) 

Combining this measurement of mode lifetime with the previ- 
ous one, and assuming that they are independent estimates, we 
finally get: 

log 10 (T/d) = 1.21^ (4) 

or 16 d (1-cr range between 9 and 50 days). Note that the dis- 
tribution about the mean is very asymmetric which implies that 
the upper limit is quite tentative. Moreover, the quoted lifetime 
is very close to the bin width (resolution) in the power spec- 
trum and the mean can thus be influenced by the observation 
time length. However, the lower limit of 9 days is quite secure 
and valuable. 

An interesting work to test the accuracy of the lifetime 
would be to compare the lifetime deduced from the scatter of 
the I — and 1 = 2 modes. Unlike t = modes, ( = 2 modes can 
be affected by the rotational splitting (vsini of 2.6 km s" 1 de- 
termined from the CORALIE spectra) and their scatter should 
be larger. However, the number and the non-consecutive po- 
sition at all of such modes (see Fig. [3) do not allow us to do 
it. 

3.3. Frequency accuracy 

Now that we have determined the mode lifetime, we can esti- 
mate the uncertainties in the mode frequencies. To do this, we 
made 850 additional simulations using several different ampli- 
tudes and assuming noise and sampling similar to the v Indi se- 
ries in order to determine the scatter of the detected modes. The 
result of the simulations allowed an estimate of the accuracy for 
the individual detected frequencies, including the effect of the 
mode lifetime, and these are given in Table [2] 



of the finite lifetime is to randomly shift each oscillation peak 
from its true position by a small amount. This is responsible for 
the scatter of the frequencies. 

The scatters are determined by assuming that peaks are un- 
affected by the rotation of the star and that the large spacing 
varies only slightly between three consecutive detected modes. 
For each measured I =0 frequency, except those at the ends of 
the ridges, we calculated the difference between the measured 
frequency and that expected from linear interpolating the po- 
sitions of the two nearest neighbors. This quantity was then 
transformed to the rms scatter of that peak about its expected 
position, in order to easily compare it with simulations (see 
Kjeldsen et al. 120051 1 . In this way, three values of scatter per 



4. Conclusions 

Our observations of v Indi from two sites have allowed us 
to identify 13 oscillation modes. The mode identification was 
complicated by the apparent effects of avoided crossings. We 
found a large and small separation of Av = 25.14 + 0.09juHz 
and 6v 02 = 2.96 + 0.22 fiHz at 330 //Hz. 

Based on the scatter of the observed frequencies, we in- 
ferred a mode lifetime of about 16 days (1-cr range between 9 
and 50 days). The lower limit of 9 days is secure whereas the 
mean and upper limits values can be influenced by the observa- 
tion time length and have to be taken with caution. Note that the 
lifetime of v Indi, even its lower limit of 9 days, is far greater 
than previous detections on other stars as the Sun (Chaplin et 
al. [T997] l, the dwarf stars a Cen A and B (Kjeldsen et al. 120051 
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and the giant £ Hya (Stello et al. 2006 ). This new measurement 
of damping time, which is until now a poorly known parame- 
ter, is very different from the previous ones and thus present a 
challenge to theoretical models. The theoretical study of v Indi, 
with asteroseismic and non-asteroseismic constraints, is post- 
poned to a third paper. 
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